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We report on the discovery of the simultaneous observation of Fishtail Effect (FE) and Peak Effect (PE) via the study of 
angular dependence of dc magnetization hysteresis loops in a clean crystal of 2H-NbSe2. These result clarify and establish the 
occurrence of the reentrant characteristics in order to disorder transformation in an isothermal scan close to zero field 
superconducting transition temperature of 2H- NbSe2. The ubiquitous FE arises from coalescence of two anomalous 
variations in current density J c , one of which is related to the collapse of the rigidity FLL close to H C 2 and another one 
located at very low fields presumably corresponds to the pinning induced order to disorder transformation. When the two 
effects get well separated, an ordered state of FLL exists in the intermediate field region. 
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The vortex state of type II superconductors conceived 
by Abrikosov [jjj encompasses the translational symme- 
try and such an ideal flux line lattice (FLL) is expected to 
move with arbitrarily small transport current under the 
influence of Lorentz force. The chemical impurities and 
structural imperfections in the underlying atomic lattice 
in any real sample provide the sources for preferentially 
pinning the flux lines at the expense of destroying the 
perfect translational symmetry of the FLL and thereby 
impart a critical current density (J c ) to the vortex ar- 
ray. J c of any superconducting specimen is expected to 
monotonically decrease while approaching the supercon- 
ducting to normal phase boundary as the field increases 
(at fixed T) or as the temperature increases (at fixed H). 
However, experimentally one encounters the phenomenon 
of anomalous maximum in J c [^]-[5| with increase in H(or 
T) in all varieties of superconducting samples. This 
anomalous maximum in J c is referred to either as Peak 
Effect (PE) @ or as Fishtail Effect (FE) @. The basis of 
the latter nomenclature namely, FE, is the characteristic 
shape of the isothermal dc magnetization hysterisis loop 
of a pinned type-II superconductor. 

As the concept of pinning relates inversely to the no- 
tion of perfect periodicity in the underlying symmetry 
of FLL, it is now generally accepted that the anoma- 
lous maximum in J c elucidates the occurrence of or- 
der to disorder transformation [|Hj]-[l2| in the FLL as a 
consequence of competition and interplay between the 
elastic energy, pinning energy and thermal fluctuations 
fl3|| . The advent of high T c superconductors focused 
widespread attention Jl3[ ] on the issue of thermal melting 
E3 of the pure Abrikosov FLL. In recent years, convinc- 



ing data have been collated in favor of first order nature 
of melting of FLL in very pure and nearly stoichiometric 
single crystals of High T c cuprates [[f5]-[l8| . In the well 
studied YBa2Cu30y system, there are evidences that a 
sharp PE, exists in juxtaposition to the FLL, melting 
transition ]l7| and this coincidence supports a widely 
held view that a sharp PE located at the edge of de- 
pinning transition of FLL is a finger print of collapse of 
the elastic moduli of the FLL. In slightly off- stoichio- 
metric single crystals of YBa2Cu307 S, the maximum 
in J c has been reported to evolve continuously from a 
narrow peak lying close to H C 2 (i.e., the usual Peak Ef- 
fect) to a broad hump extending over a large field region 
located far away from H C 2 (i.e., the Fishtail Effect) as 
the effective pinning increases as a consequence of ei- 
ther the increase in quenched random disorder or the 
decrease in temperature. Thus, notionally, the PE and 
the FE may be treated as distinct and mutually exclu- 
sive features in dc magnetization hysterisis data though 
intimately related. In the context of single crystals of 
Bismuth Cuprate system, in which there have been re- 
ports of only the fishtail type of anomalous maximum in 
J c jfJI, it was demonstrated from /iSK studies that 
the spatial order of FLL undergoes a sharp change at a 
field value coincident with the onset of fishtail anomaly 
in dc magnetization hysteresis data. Another fiSK PQ] 
study in the same Bismuth Cuprate system in different 
(H, T) region had established that the spatial order of 
the FLL undergoes a sudden reduction at (H, T) val- 
ues where melting of FLL was anticipated. Thus, both 
peak effect and fishtail effect could be associated with 
the notion of occurrence of order to disorder transforma- 
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tion which is either first order or second order. There 
have been numerous reports in high T c superconductors 
wherein different phase boundaries across which vortex 
state undergoes order to disorder transformations corre- 
sponding to onset and peak position of fishtail effect have 
been drawn in the (H, T) space. Some of these phase 
boundaries match with those phase boundaries drawn 
from other transport and thermodynamic measurements 
whereas some others appear to approach them at multi- 
critical points (see, for instance, Fig. 4 in (|] and Fig. 4 in 

Considering that different phases of vortex matter and 
transformation amongst them are a consequence of com- 
petition and interplay between elastic energy , pinning 
energy and thermal fluctuations, it is of interest to ex- 
plore the presence of such characteristics in the vortex 
states of conventional low temperature superconductors 
in the appropriate (H,T) region. Amongst the low tem- 
perature superconductors, the clean single crystals of 2H- 
NbSe2 system having T c (0) « 7 K arc in current focus 
pjj , from the point of view of pristine physics issues of 
vortex matter |L3| because of (i) their very weak pinning 
nature, (ii) their relatively large Ginzburg number value 
(Gi ~ 10 ^) J@L an d (iii) the existence of a robust PE 
in them [p|, p2p3| . We present here new results pertain- 
ing to the discovery of simultaneous presence of fishtail 
effect and peak effect in a clean crystal of 2H-NbSe2 at 
low fields (in H < 1 kOe where the FLL lattice constant 
a > 2000 A°) and in the high temperature (T/T c (0) > 
0.95) region. We believe that these results are the first 
of their kind. They pertain to a vortex array in its di- 
lute limit where the interaction amongst the vortices is 
in a nascent stage and the fluctuation effects are very 
strong due to the close proximity to T c (0). These results 
have lead us to construct phase boundaries in (H,T) space 
which have the potential to clarify complex issues emerg- 
ing from studies related to PE/FE/FLL melting in high 
T c cuprates p^-|l8|. One specific issue that the simulta- 
neous observation of FE and PE aims to resolve is the 
reentrant characteristic of order-disorder phase boundary 
of FLL as the interaction gets progressively 

enhanced (with increasing field) while thermal fluctua- 
tions (T) and quenched random disorder remain fixed. 

For our present study we chose a single crystal of 2H- 
NbSe 2 (T C (0)~7.2K with a AT C ~50 mK). The sample 
dimensions are 5X2X1 mm 3 and J c values in it lie be- 
low 1000 A/cm 2 . It is our belief that the present sample 
has a level of purity which lies in between those of the 
crystals A and B, used by us in our earlier studies [^) . We 
measured its ac magnetization screening response using 
a high sensitivity ac-susceptometcr p5fl and dc- magne- 
tization using a standard Quantum Design SQUID mag- 
netometer with a specially designed home made sample 
holder which allows an angular variation of 0° to 180° 
between the field H and the ab-plane of the single crystal 

Fig.l shows the temperature dependence of the dia- 
magnetic screening response (x'(T)) for the FLL created 
at different (|| c) values. From a generalized Critical 



State Model f|| relation [§7) : 



, where a is a geometry and size dependent factor and 
J c is the (H,T)- dependent critical current density. From 
this relationship, it can be deduced that PE phenomenon 
(i.e., peak in J c ) should manifest itself as an anomalous 
increase in the diamagnetic screening response (x'(T)) 
for a given Ebj c . In Fig. 1, we have marked with arrows 
the peak temperatures T p of the PE in x'(T) behavior at 
various H<fc. It can be seen that T p values increase with 
decreasing H<2 C down to 200 Oe. Further more, another 
noteworthy feature in x'(T) data is that, the PE starts 
to broaden substantially about T p value as one moves 
to lower fields (from 100 Oe to 50 Oe). To reveal the 
possible connection between the phenomenon of broad- 
ening of PE and the behavior of FLL at low fields it may 
be pertinent to point out here that, in the field range 
from 100 Oe to 50 Oe where the PE starts to broaden, 
the values of a vary from 4800 A° to about 6800 A°. 
Recalling that A c ~ 4550 A°, in 2H-NbSe 2 in the (H,T) 
region under consideration |2{| , the vortex array is in the 
dilute limit (ao > A) for H < 100 Oe and the flux lines 
are only weakly interacting , with the result the FLL is 
in a state which is easily susceptible to thermal fluctua- 
tions and pinning effects. The inset in Fig.l shows the 
PE curve which has been determined by picking out the 
field - temperature values at which the peak in PE oc- 
curs in x'(T) data. The noteworthy features of PE are : 
(i) at high fields (H > 200 Oe) the PE curve tracks the 
H C 2(T) curve, such that both of them have a slope of ~ 
5 kOe/K, however at lower fields (H < 100 Oe) the PE 
curve bends away from the H C 2(T) curve, (ii) the pro- 
gressive increase in the size of the error bars on the data 
points indicates the commencement of process of broad- 
ening out of PE phenomenon, while approaching T c (0). 
Our previous studies || had revealed that increased ef- 
fective pinning broadened out the PE, we thus surmise 
that an interplay between pinning and thermal fluctua- 
tion effects at low H - high T determines the broaden- 
ing feature in PE in Fig.l. This prompts to investigate 
the features of PE and its broadening through isother- 
mal magnetization M(H) hysterisis loop measurements, 
where we study the progression of flux line lattice of vary- 
ing ao under the influence of fixed thermal energy. To 
perform these isothermal M(H) measurements, we chose 
particularly those temperatures (close to T c (0)), where 
the PE broadened and the PE curve seems to bend away 
from H c2 (T) (see Fig.l). 

A superconductor shields itself from any change in 
external magnetic field by setting up currents equal to 
J C (B), which in turn result in the magnetization M at 
a given field(H). When one obtains a dc magnetization 
hysterisis loop, the width of the magnetization hysteri- 
sis at a given field, i.e., AM (AM(H)= M(H1)-M(H|)) 
may be taken as a measure of J c . This implies that any 
non- monotonic variation in the behavior of J C (H) can, 
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therefore, show up as an anomalous modulation in the 
width (AM) of the magnetization hysterisis loop. The 
PE emerges thus as an anomalous bubble like anomaly 
superposed on the quasi reversible magnetization (M(H)) 
hysteresis loop (see PE region centered around H p in the 
main panel of Fig. 2(a)). In isothermal dc magnetization 
hysterisis data at T=6.95 K in 2H-NbSe2, we observe 
that there is a well formed PE (anomalous M(H) hyster- 
isis bubble) centered around H p «1000 Oe (cf. Fig. 2(a)) 
and this peak field value is consistent with the T p (H) 
curve in Fig.l. In Fig. 2(a) in the field range from 50 Oe 
to 250 Oe (see arrow marked in Fig. 2 (a)), there is 
an anomalous modulation in the width of the magnetiza- 
tion hysterisis loop. We can determine the J c using the 
relationship J c cx AM/d, where d is the thickness of the 
sample. In the inset of Fig. 2(a) the behavior of J c versus 
H at T=6.95 K is shown in dark square symbols on a 
semi- log plot. Such a J C (H) plot shows an overall linear 
(i.e., J c oc exp(— H)) behavior (examine the straight line 
passing through the data points). It can be easily viewed 
that in the field range of 70 Oe < H < 250 Oe, there is an 
anomalous behavior in J C (H) which brings out the devia- 
tion from the overall linear behavior. We label the center 
of gravity of this anomaly in J c as at » 100 Oe. The 
development of this anomaly at H^ as a function of tem- 
perature is investigated by studying the M(H) hysterisis 
loop (cf. main panel of Fig. 2(c)) at T = 7.0 K, where the 
hysteresis loop seems to be completely anomalous and 
irreversible in the entire field range from Oe < H < 
H c2 . The visible differences that emerge by comparing 
the M(H) hysterisis loop at T = 7.0 K and T = 6.95 K 
are : (i) one cannot precisely locate the low field anomaly 
at Hd in the M(H) loop at 7.0 K and (ii) one cannot also 
clearly identify from the M(H) hysteresis curve the upper 
PE which occurs close to H C 2 . In the inset of Fig. 2(c), 
we have put together the behavior of the J C (H) at 7.0 K 
and 6.95 K. The usual expectation is J C (H, T = 7. OK ) 
< J C (H, T = 6.95 K), however such an inequality is sat- 
isfied only for H < H d . At H > H d , the J C (H, T = 7.0K 
) > J C (H, T = 6.95 K ), which is a manifestation of an 
anomalous behavior in J c and it signals the occurrence 
of peak effect like feature in J c at low fields which are far 
from H C 2. It thus seems that the anomalous modulation 
in J c at at 6.95 K in Fig. 2(a) is a fingerprint of a low 
field PE phenomenon and it presumably survives at 7.0 
K. 

We now present an analysis of the pinning force 
F p (=J c *H), which brings out features not readily appar- 
ent from the M(H) hysteresis data or from the J C (H) 
behavior. We first determine F p using the J C (H) data 
extracted from the M(H) hysterisis loops, and we call it 
F p (measured) . Next, we construct a monotonically de- 
creasing J C (H) which is devoid of any anomalous mod- 
ulation, i.e., this constructed J c gets rid of any PE like 
features at low fields or at high fields (a procedure to con- 
struct such a J C (H) can be seen from the semi log plot 
of J C (H) (in the inset of Fig.2(a)) at 6.95 K). Using this 
constructed J c , we determine the pinning force which we 



call F p (constructed) =J C (constructed)* H. Finally, we de- 
termine AF p =F p (measured) — F p (constructed). AF V is 
a quantity which can convey signatures of anomalous |29| ] 
changes in the intrinsic pinning force density being expe- 
rienced by the flux lines. In Figs. 2(b), 2(d), 2(f), 2(h), 
we have plotted AF p as a function of H. At 6.95 K (cf. 
Fig. 2(b)), one can clearly distinguish two well resolved 
PE features in AF p (H). One PE located at H d -100 Oe, 
is the novel low field PE, this peak in AF p at H^ coincides 
with the field at which the anomalous non monotonic low 
field behavior in J C (H) is observed (cf. inset of Fig. 2(a)). 
The other PE, is the high field PE at H p , which occurs 
close to H C 2. At T=7.0 K (cf. Fig. 2(d)), one can again 
distinguish two well resolved peak like features in AF p , 
at ELi and H p . It should be noted that although the M(H) 
hysteresis loop at 7. OK was devoid of any distinguishing 
features of PE at H^ and H p , through the present anal- 
ysis we can now discern features of PE at both H^ and 
H p in the AF p (H) plot at 7.0 K. If one now compares the 
Fig.2(b) with Fig.2(d), it seems that though the PE at 
Hrf does not move as we change T from 6.95 K to 7.0 K, 
the H p has moved from -1000 Oe at 6.95 K to ~750 Oe 
at 7.0 K, which is a rate of 5 kOe/K. This slope value of 
dH p /dT, matches well with the slope of the PE curve in 
the inset of Fig.l, which was determined from temper- 
ature dependent ac susceptibility measurements. From 
the above analysis one can conjecture that due to a mo- 
tion to lower fields (close to H^), of the PE at H p at the 
rate of 5 kOe/K as T is increased, one should be able to 
see the emergence of a single peak from the coalescence 
of two well resolved PE peaks in AF p (i.e., one at H^ and 
another at H p ). We now show the M(H) hysteresis loop 
at subsequently higher temperatures. The M(H) hystere- 
sis loop at T=7.05 K (cf. Fig. 2(e)), appears even more 
anomalous than that at 7.0 K. It is a broad, irreversible 
loop which results in a single broad peak in AF p (cf. Fig. 
2(f)). The M(H) hysterisis loop at 7.05 K now re- 
sembles so called "Fishtail Effect Anomaly". This 
fishtail effect like anomaly in the M(H) hysteresis loop 
and in AF p (H) plot persists at T > 7.05 K. Figs. 2(g) 
and 2(h) show M(H) hysterisis loop and the F p (H) at T = 
7.1 K. It is to be noted that the width of the FE anomaly, 
(i.e., AH d in Fig.2(b), Fig.2(d), Fig.2(f), Fig.2(h)) has 
an interesting temperature dependence, we shall discuss 
it at a later stage. We shall now demonstrate through 
a sequence of M(H) hysterisis loops obtained at various 
orientations with respect to the ab plane of the single 
crystal as to how the composite broad FE effect at 7.0 K 
(see Fig. 2(c)) gets resolved into two anomalous maxima 
in J c which are centered around H^ and H p . 

It had been argued by Pippard ||(J and also ex- 
perimentally seen that the PE peak H p always scales 
with H C 2. 2H-NbSe2 is an anisotropic system to 
which anisotropic Ginzburg-Landau description ap- 
plies pi} ). The angular dependence of H c2 is given 
by as H c2 =H C 2(||c,T)(Sin 2 (6')+e 2 Cos 2 (e'))- 1 / 2 , where 
e=l/(Anisotropy of crystal) and 9 the angle between H 
and the ab plane of the single crystal of NbSe2. As 6 
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changes from 90° to 0°, H C 2 increases from the value of 
H C 2(||c,T) to that of H C 2([|ab,T), and concomitantly the 
value of Hp also follow suit. In Fig. 3, the M(H) hysterisis 
loops in 2H-NbSe2 at 7.0 K for different The insets in 
each of the Figs. 7(a) to (f) show F p vs H/H c2 (0). In the 
inset of Fig. 3 (a) we can see that the low field, maximum 
in F p at Hd and the high field PE peak at Hp lie in jux- 
taposition. When the angle 9 changes from 90° towards 
0°, the following features in Fig. 3 are noteworthy: (i) In 
all the insets of Fig. 3, the peak in F p at H p seems to oc- 
cur at H/H c2 (0)«O.7. (ii) The low field peak in F p at H d 
shifts to lower values of H/H C 2(0) as increases (iii) The 
most striking feature emerging from the plots of Fig. 3 is 
the separating out of the two anomalous modulations in 
J c corresponding to a distinct double peak structure in 
the Fp(H) plot. The main panels of Fig. 3 show that by 
varying 0, the PE peak at H p moves to higher fields (since 
Hp/H c2 (0)«O.7), and as the H p value moves sufficiently 
far away, the novel low field anomalous modulation in 
J c at Hd gets identified distinctly as an independent fea- 
ture. If we attempt to fit the observed H c2 (0) behavior to 
the above stated Ginzburg-Landau relationship, we get 
a value of e as ~0.5, or an anisotropy of about 2 for 2H- 
NbSe 2 , which is a satisfactory value for NbSe2 for H < 2 
kOe ps| . Thus our anisotropy study of the M(H) hyster- 
isis loop elucidates the fact that, FE could be composed 
of an anomalous modulation in J c centered at low field 
Hd and the usual PE peak at H p . 

We now summarize our discussion and propose a 
schematic for the low H and high T part of the vortex 
phase diagram. In the inset of Fig. 1(a), we had presented 
a phase diagram which comprised the PE phase bound- 
ary T p (H) across which the FLL disorders into an amor- 
phous phase. However, this phase diagram is incomplete 
as it does not give us any indication of changes in the 
vortex matter occurring at field values such as at H^ and 
AHd, which we can observe from the isothermal M(H) 
hysteresis loops. Returning to Fig. 2, we consider the be- 
havior of the width AH^ of the low field PE centered at 
Hrf. Shown in Fig. 4, is a shaded region which is obtained 
by estimating the full width at half maximum of the peak 
in AFp at H d in AF p (H) plots (cf. Fig.2) at different T's. 
This high pinning shaded region is bounded by two curves 
R d lower and H/w er , where AH d = \R d u PP er - U d lower \. 
The width AH d seems to increase with T and reaches 
a maximum at about T = 7.05 K and above this tem- 
perature (cf. Fig.2), the width AH^ appears to decrease 
again. At very low H (ao > A), where the rigidity of 
the FLL (viz. C66 ~ exp(-ao/A)) is exponentially small, 
the FLL is highly susceptible to thermal fluctuations and 
pinning effects which could disorder the FLL. For vortex 
states corresponding to the shaded region of Fig. 4 (low 
H - high T), the weak rigidity permits thermal fluctu- 
ations and pinning effects to dominate as evidenced by 
the broad hump in the M-H curves and hence the FLL 
states in the shaded region belong to "disordered phase" 
of vortex matter which is presumably characterized by 
a small Larkin volume V c . We now attempt to explore 



features in our temperature dependent ac susceptibility 
and dc magnetization data which can possibly be linked 
to the phase labelled as the "disordered phase" in the 
H-T plane . As implied from the above analysis, AH^ 
identifies the field span over which J C (H) varies anoma- 
lously. In temperature dependent ac susceptibility data, 
the broadening of PE peak commenced in isofield runs 
at H < IkOe. We believe that the said broadening is a 
consequence of encountering the shaded region in fixed 
H scan (see Fig. 4). On the other hand, in fixed T scan 
in dc magnetisation hysteresis data, the "lower" PE and 
"upper" PE appear as distinct features if H d u PP er line lies 
sufficiently far below the H C 2 line. In the field range be- 
tween the H £ j ttpper and the PE region, the vortex state do 
not show any anomalous variation in J c . We, therefore, 
feel that such FLL states identify the well ordered vor- 
tex solid. However in the field region above the point A 
(marked in main panel of Fig. 4 ), where the shaded region 
and the H p line have merged, if one records an isother- 
mal M(H) hysteresis loop, then this would be the region 
where the two PE would have merged and would have 
the appearance of the Fishtail like effect. Above H p and 
below Hi rr curves || , the FLL is presumably in an amor- 
phous state with non zero pinning. Such amorphous state 
maybe indistinguishable in its pinning behavior from the 
"disordered phase" in the shaded region. Hence, we pro- 
pose that the vortex region between H p and Hi rr lines 
and the shaded region labeled as the "disordered phase" 
essentially identify the same vortex phase. In the inset 
of Fig. 4, we show the phase boundary across which the 
ordered solid disorders. This has been constructed by 
combining the H p curve, to the Yi d upper curve (see main 
panel of Fig. 4). This new phase boundary across, desig- 
nated as Hord in the inset of Fig. 4, has a reentrant nature, 
i.e., in fixed T scan, the order-disorder boundary shall 
be encountered twice. Such a reentrant characteristic is 
reminiscent of reentrant nature of the FLL melting curve 
of a pure (pinnin g fr ee) Abrikosov state, first proposed 
by D. R. Nelson [[13|. Between this newly constructed 
reentrant H ord (T) line and the Hi rr (T) line, the vortex 
matter can be designated as the "pinned liquid phase" , 
following the nomenclature in vogue in the literature on 
vortex state studies Jp| . 

In conclusion we may state that 2H-NbSe2, seems to be 
an appropriate system, whose intrinsic parameters, make 
the phase boundaries, across which the various phases 
of vortex matter exists, become accessible to within ex- 
perimental limits. Though the present set of results are 
the first preliminary steps in characterizing the different 
phases of vortex matter, the experimental and theoreti- 
cal investigations in this field are wide open to enormous 
possibilities and opportunities. 

We would like to acknowledge Nandini Trivedi, Satya 
Majumder and Mahesh Chandran for fruitful discussions. 
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FIG. 1. Temperature variation of ac susceptibility (f=211 
Hz, h oc =l Oe(r.m.s.)) in a single crystal of 2H-NbSe2 for Hd c 
|| c. The peak temperature T p have been identified by arrow. 
Note that PE peak broadens as H^c progressively decreases 
for Hdc < 1 kOe. The inset shows the plot of T P (H) and 
T C (H) values in the low field region of the magnetic phase 
diagram. 

FIG. 2. Isothermal dc magnetisation loops, extracting J c 
values and computed pinning force values in 2H-NbSe2 for 
Hdc||c at various temperatures. For explanation of various 
symbols used in this figure, see text. 

FIG. 3. Angular dependence of magnetic hysteresis data 
in a crystal of 2H-NbSe2 at 7.0 K. The angle 9 is measured 
w.r.t. basal plane of hexagonal atomic lattice of 2H-NbSe2. 
For explanation of various symbols in this figure, see text. 

FIG. 4. Magnetic phase diagram corresponding to tem- 
perature dependent ac susceptibility (Fig.l) and field de- 
pendent dc magnetisation hysterisis (Fig.l) data. The inset 
shows, schematically the proposed magnetic phase diagram 
constructed the results of ac and dc magnetisation experi- 
ments. See text for explanation of different phases bound- 
aries. 
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